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The single particle spectral-weight function (SWF) of the ionic Hubbard model 
at half filling is calculated in the cluster perturbation theory approximation. An 
abrupt change of regime in the low-energy region, near the chemical potential, is 
found at a critical value, Uc, of the coupling constant (Hubbard U). The SWF at 
the Fermi points kF=±TT/2 jumps, as U increases, from a two-peak structure, the 
gap edges, to a four-peak structure accompanied by a (non-vanishing) minimum of 
the charge-gap. The two inner peaks of this structure show very small dispersion 
(flat bands) away from the Fermi points, whereas the outer peaks mark the edges of 
the Hubbard bands. No other signatures of abrupt change are detected in the SWF. 
The two regimes are physically realized in the angle-resolved photoelectron spectra 
of {TaSe4)2l , and the blue-bronze Kq^^MoO-^, respectively. 

PACS numbers: 71.10.-w, 71.10.Fd, 71.10.Hf, 71.30.+h 



I. INTRODUCTION 



Quasi-one-dimensional (QID) systems have been the object of intense experimental and 
theoretical activity over the last twenty years. They show highly anisotropic properties, 
with a privileged direction of enhanced charge transport. Their interest lies in the hope that 
they can be good candidates for the physical realization of non-Fermi liquid behavior. This 
interest, in low-D systems in general, has expanded very rapidly in recent years partly due 



2 



to the technological development of low-D artificial structures and nano-scale materials. 

Above their Peirls temperature (or when doped away from half filling), these QID sys- 
tems are conductors and display Luttinger liquid behavioi>i, i.e., the absence of quasi-particle 
excitations in the Fermi liquid sense (a quasi-particle peak at the Fermi level) and the 
excitation, instead, of decoupled collective modes of charge (holons) and spin (spinous) 
character, a phenomenon usually known as charge-spin separation. The absence of a 
Fermi edge has indeed been found in angle-resolved photoelectron spectroscopy (ARPES) 
of {TaSe4)2l^'^, Kq^^MoO^'^, and the organic conductor TTF-TCNQ (tetrathiafulvalene- 
tetracyanoquinodimethane)^'^. For a good review, see Ref 7. Clear experimental signatures 
of spin-charge separation are, however, very scarce in QID conductors, with the notable 
exception of TTF-TCNQ reported very recently^. All these compounds in their metal- 
lic state can be modeled by the low-energy physics of the doped ID single-band Hubbard 
Hamiltonian^. Alternatively, they have been analyzed on the basis of the Luttinger model or 
the Luther- Emery model (when a spin gap is expected) . Some puzzles still remain unsolved^ 

Below their Peirls temperature, these QID compounds, as well as many others like 
halogen-bridged transition-metal chains, conjugated polymers, and organic charge-transfer 
salts are usually insulating. In these systems, the competition between strong on-site cor- 
relations and their kinetic energy gives rise to significant localization of their itinerant elec- 
trons. This often leads to the stabilization of non-metallic ground states with or without 
charge- density waves (CDW). Thus {TaSe4^)2l and the blue bronze Ko^MoOs, for instance, 
are CDW insulators^, while the nearly ideal ID CuO chains in SrCu02 and Sr2CuO^ are 
responsible for the insulating character of these charge-transfer insulators^. Signatures of 
spin-charge separation have also been found in ARPES of SrCu02^ and in the dielectric 
response of Sr2CuO^. 

Most of these insulating systems can be conveniently described by the Emery modeiii, 
which is a generalized Hubbard Hamiltonian on a two-sublattice model, made of cations 
(say Cu ions) and anions (say O ions), respectively. The on-site energy levels and repulsions 
{^d,Udd) and {ep,Upp), for Cu d orbitals and O p orbitals, are coupled by nearest-neighbor 
hopping of strength tdp and Coulomb repulsions Vdp. Its ID version reads 



H = ndis + Udd ^ ndi^Udii+e. 
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+tdp {dtsPjs + he) + Vdp TLdiUpj 

<ij>s <ij> 

where creates an electron (hole) in a ci orbital at site i with spin s, i running over all 
the Cu sites. Similarly pjs anihilates an electron (hole) in a p orbital at site j and spin s, 
j running over all the O sites. As usual, < ij > means summation over nearest-neighbors, 
fig = c'^Cs and = + rii^ denotes the charge at the ith site. 

We may now simplify this hamiltonian. At first sight the Cu-0 repulsion seems essential 
since a large charge-transfer is expected. However, Vdp is usually much smaller than the 
on-site repulsions and, furthermore, this term gives rise to new physics only in the event of 
exciton formation. Hence, if we are not especially interested in these processes, V^p can be 
safely ignored. One is then left with a charge-transfer model Hamiltonian which can still 
give a reasonable description of some of these compound a^*^ ', e.g., SrCu02 and Sr2CuO^. 
If we are now willing to reduce the number of parameters by putting (somewhat arbitrarily) 
Upp = Udd, then the so-called ionic Hubbard model (IHM) follows. It can be written simply 
as 

H = -tJ2 4c,. + y J^i-iyn^s + UY1 ^^T^^i (2) 

<ij>s is i 

where A is the on-site energy difference between even and odd sites, usually known as the 
charge-transfer energy. 

Although, strictly speaking, this Hamiltonian does not describe accurately any specific 
system (since quite generally Upp < Udd), it provides a simple, minimal, model where the 
interplay among covalency (t), ionicity (A) and correlation (U) gives rise to a rich phase 
diagram within which different ID compounds can be placed. Originally proposed by Na- 
gaosa and Takimotc^^ as a model for ferroelectric perovskites and later by Egami et al^^*^ 
to explain the neutral-ionic transition in some organic crystals, this Hamiltonian is ideal 
for studying the nature of quantum phase transitions in ID electron systems. On general 
grounds, one expects a transition from an ionic, weakly-correlated band insulator (BI) phase 
to a neutral, strongly-correlated Mott insulator (MI) phase as U increases. An important 
and controverted issue is the nature of this transition as well as whether two critical points 
rather than one separate both phases. Depending on the method of calculation used, either 
on o^7ii?ii? or two22*21i22i2i critical points have been predicted, so that the controversy cannot 
be considered as closed yet. 

In this paper we stay outside this controversy and will rather concentrate on the single- 



4 



A(k,E) 




I 1 1 1 1 

-2-1012 

E-n(eV) 

FIG. 1: Single-particle spectral-weight function A(k,E) for the half-filled ionic Hubbard model at 
the Fermi points k = ib7r/2. From top to bottom, [7 = 1, 2.0775 (just below Uc), 2.078 (just above 
Uc), and 3. We have taken t = —1 and A = 1. All the energies are given in eV. 

particle spectral-weight function (SWF) A{k,uj) of the ID IHM at half filling, which can be 
compared with ARPES of several insulating materials. A{k,uj) will be calculated using the 
cluster perturbation theory (CPT) approach of Senechal et aF'*. This is briefly described in 
Sec II. Sec III gives our results for A{k,u!) (a sudden change of regime at a critical value of 
U) followed by a discussion showing that {TaSei)2l and f^o.sMoOs, are good examples of 
these two regimes. Finally, Sec IV closes the paper with some concluding remarks. 

II. CLUSTER PERTURBATION THEORY (CPT) 

Since this method has been discussed at length in Ref 24, we simply summarize it very 
briefly here. In CPT one divides the lattice (here the chain) into a number of equal clusters. 
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The single-particle Green's function (GF) on these clusters is then found by exact diago- 
nalization with open boundary conditions. We have made use of a variant of the Lanczos 
algorithm specifically designed to calculate dynamic quantities-^. The approximation now 
consists in neglecting the intercluster self-energy, so that the GF's of neighboring clusters 
are connected by hopping terms only. Periodic boundary conditions are then imposed on 
the whole chain, i.e., between the extreme clusters. To be specific, let mi denote the site m 
of the cluster i. The exact Greens function Gmi,nji of the whole chain is given by the well- 
known Dyson's equation in matrix form (Gg ^ — S)^ = /, in terms of the non-interacting GF 
and the exact S. In CPT this exact S is approximated by Smi.nj = where is the 

cluster self-energy matrix. This approximation is applicable to any lattice in any dimension. 
It can be understood as a lowest-order contribution to a systematic perturbation expansion 
in powers of the intercluster hopping^^^. It turns out, on the other hand, that CPT is a 
limiting case of a more general variational cluster approacb^I. 

In this paper we concentrate on the SWF A{k,uj), given as usual by 

A{k,uo) = --\mG{k,uj + 11^) (3) 
7r 

where G{k, uj) is the Fourier transform (FT) of the single-particle retarded GF and r] a small 
positive number. This FT must be calculated with some care since Gmnii — j) is periodic 
in ij but not in mn due to the open boundary conditions used in the clusters. The correct 
formula ia^^^ 

G{k, ^) = ^ E e~^'('"-"^G„n(iVA;, u) (4) 

mn 

where N is the number of sites in a cluster and Gmn{k) the FT of GmrM ~ j)- 

III. THE SPECTRAL WEIGHT FUNCTION OF THE IONIC HUBBARD 

MODEL AT HALF FILLING 



We consider a chain of ninety-six sites of two kinds with levels at ±A/2, at half filling, 
and take t = —leV and A = leV. Clusters of eight sites have been adopted after checking 
that increasing the cluster size up to twelve sites does not change much the results for the 
whole chain. Fig 1 shows A{k,E) at the Fermi points kp = ±7r/2 for increasing U. A 
broadening r] = 0.05 eV has been given to the otherwise delta functions. Two regimes are 
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FIG. 2: Density of states for the same model of Figl with U = 2.0775, eV just below Uc- 

observed in the low-energy region: For small U, a two-peak structure is seen close to the 
chemical potential fi (zero reference-energy). This structure persists, while the two peaks 
approach each other, up to t/ = 2.0775eV. Abruptly, at U = Uc = 2.078 eV a four-peak 
structure appears, with the peak heights being roughly half those of the previous two-peak 
structure. It seems, therefore, that each of the latter peaks splits, developing a new peak 
at higher energy and thus leading to the new four-peak structure. As U further increases, 
the inner peaks, closer to the chemical potential, start to separate from each other quickly 
approaching the outer peaks which in turn move apart more slowly. Finally, for very large 
U, much larger than t and A, the SWF tends asymptotically to a Mott-Hubbard situation 
without any signature of abrupt change. The single-particle (charge) gap, after passing 
through a minimum at Uc, increases again. It never vanishes, in agreement with recent 
findingsi^i^i. 

A somewhat different, but complementary perspective is afforded by the density of states 
(DOS), N{E) = (1/M) J2kMk,E), where M is the number of k's. Fig 2 shows the DOS 
for U = 2.0775eV^, just below Uc- It is reminiscent of the non-interacting DOS (a brunch- 
cut surrounded by two singularities and cut by a gap at the chemical potential). This 
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FIG. 3: The same as Fig 2 but for U = 2.0780, eV just above Uc 



region U < Uc, is the band insulating regime, the DOS being similar to that for U — 0, 
but with a continuously decreasing gap. Fig 3, on the other hand, displays the DOS for 
U — 2.078ey, just above in correspondence with the third panel (from top to bottom) 
of Fig 1. The two outer peaks in this latter figure mark the top (bottom) of the occupied 
(empty) Hubbard band. The two inner peaks stand isolated near the chemical potential 
resembling two "impurity-like", non- dispersive peaks within the (wider) gap between the 
Hubbard bands. All through the region U > Uc, the DOS is quite different from that at 
C/ = 0. 

To check what changes have taken place at A; 7^ ±7r/2. Figs 4 and 5 display A(k,E) for the 
same U values as the above DOS along the large Brillouin Zone (BZ), — tt < A; < tt (in the 

extended zone scheme). The portion | k \> 7r/2 can be mapped, if one wishes, onto the small 
BZ giving the empty bands for | k \< 7r/2 above the occupied bands. An offset has been 



8 




E-n(eV) 



k=7i/2 



k=0 



k=-7i/2 



k=-7i 



FIG. 4: Spectral-weight function A{k, E) for the same situation of Fig 2. An offset has been 
given to the plots for different k's along the large Brillouin zone (— vr < A; < vr) in order to avoid 
superposition. The figures along the left vertical axis count the number of k's starting from — tt. 
On the right vertical axis, some especial k's are indicated 



provided to the different plots to avoid superposition. The figures along the vertical axis 
number the k's starting from k = — tt (only a selected set of thirty-two k's have been shown 
for clarity). Fig 4 shows a cosine- like band cut in two by a gap at the Fermi points (cfr. 
with the second panel in Fig 1). Two shadow bands covering only part of the BZ are clearly 
visible around A; = (occupied) and k — n (empty) . They give rise to the side-peaks to the 
left and right of the main body of the DOS of Fig 2. Fig 5 should now be confronted with 
the DOS of Fig 3. The cosine-like band is now cut by a wider gap delimited by the outer 
peaks of Fig 1, third panel. The inner peaks are continued by two almost non-dispersive, 
flat bands covering only a small portion of the BZ around kp are clearly visible. The shadow 
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FIG. 5: Same as Fig 4, but for U = 2.0780, eV just above Uc 

bands around k — and k — tt have now almost disappeared along with the side-peaks, 
rather weak in Fig 3. 

It is difficult to extract the character of the peaks in a calculation based on a Lanczos 
algorithm, and much more so to trace the origin of the change of regime just described. 
However, these regimes must have some physical reality since we observe that both are 
found in ARPES of different QID insulating materials. Since most of these materials have 
occupied bands with roughly leV of bandwidth, we have accordingly taken t — —0.5eV and 
A = leV. The critical U separating both regimes is now Uc — 2eV. Let us take, just for 
illustrative purpose, {TaSe^2l, and Kq 2,MoO^, Fig 6 shows A{k, E) between k — Q and 
7r/2 for U = Q.beV < Uc, a case of the first regime, with a single occupied peak at = 7r/2. 
This one-band structure resembles that of {TaSe4)2l in the direction parallel to the chain 
(compare with Fig 2a of Ref 3). Likewise. Fig 7 shows the same information as Fig 6, but 
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FIG. 6: Spectral-weight function for (in eV) t = —0.5, A = 1 and U = 0.5 



for U = 2.1eV > Uc- We now see a two-band structure in the neighborhood of k = tt/2. 
This resembles the band structure of Kq^^MoOs (compare with Fig 15 of Ref 4). 



IV. CONCLUSIONS 

Using the cluster perturbation theory approach of Senechal et al2^^, we have calculated 
the single-particle spectral- weight function A{ku!) of the ionic ID Hubbard model at half 
filling. A change of regime is found at a critical value of U, Uc{t,A), which depends on 
both the hopping amplitude t and the on-site energy-difference, A, between even and odd 
sites. As U increases, A{kuj) jumps from a two-peak structure to a four-peak one at the 
Fermi points kf = ±7r/2. As one moves away from kp, one finds two semiconducting bands, 
for U < Uc, separated by a gap which decreases from its initial value A (at f/ = 0) down 
to a small, but non-vanishing, value at Uc- This gap is always delimited by the two-peak 
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FIG. 7: The same as Fig 6, but for U = 2.1 eV 



structure at kp- For U > Uc, instead, two flat, almost-non-dispersive bands appear around 
the Fermi points as the continuation of the inner peaks in the four-peak structure. They 
fade away very soon. The gap between them now increases from its minimum value at Uc- 
The outer peaks mark the top (bottom) of the lower (upper) Hubbard-hke bands. As U 
increases further, the flat bands approach the Hubbard bands, finally merging into them. 
Asymptotically, a Mott- Hubbard situation is approached in a continuous way, without any 
signature of abrupt change in A{k, uj) in this region of large U. 

Alternatively we can say that, for U < Uc, the two-band structure is continuously con- 
nected to the band insulator shape at U — 0. This is the band insulator regime. For U > Uc, 
two new flat bands appear and push the two wider bands of the first regime further apart. 
This is the second regime, more Hubbard-hke, which goes to the Mott-Hubbard regime at 
large U. Different QID materials have band structures which can be classified as lying in 
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either the first or the second regime. 
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